Edge electron states for quasi-one-dimensional organic conductors in the 
magnetic-field-induced spin- density- wave phases 
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We develop a microscopic picture of the electron states localized at the edges perpendicular to the 
chains in the Bechgaard salts in the quantum Hall regime. In a magnetic-field-induced spin-density- 
wave state (FISDW) characterized by an integer N, there exist N branches of chiral gapless edge 
excitations. Localization length is much longer and velocity much lower for these states than for the 
edge states parallel to the chains. We calculate the contribution of these states to the specific heat 
and propose a time-of-flight experiment to probe the propagating edge modes directly. 

PACS numbers: 74.70Kn, 75.30Fv, 73.40Hm, 73.20-r 
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Quasi-one-dimensional (Q1D) organic conductors of 
the (TMTSF) 2 X family fjj (also known as the Bechgaard 
salts) are highly anisotropic crystals that consist of par- 
allel conducting chains. The electron transfer integrals 
along the chains (in the a direction) and transverse to 
the chains (in the b and c directions) are t a = 250 meV, 
<b = 25 meV, and t c = 1.5 mcV Because of the strong 
Q1D anisotropy, the Fermi surfaces of these materials are 
open and consist of two disconnected sheets located near 
ikp, which are the Fermi momenta along the chains. In 
the presence of a strong magnetic field along the c axis, 
the interplay between the nesting property of the open 
Fermi surface and the quantization of electron orbits due 
to the magnetic field leads to a cascade of the magnetic- 
field-induced spin-density-wave (FISDW) phase transi- 
tions |^,||. FISDW opens an energy gap A at in the elec- 
tron spectrum at the Fermi level. Within each FISDW 
phase, the Hall conductivity per a-b layer per electron 
spin is quantized at zero temperature as a xy = Ne 2 /h, 
where N is an integer that characterizes the FISDW, e 
is the electron charge, and h is the Planck constant. 

The theory of the quantum Hall effect (QHE) in the 
Bechgaard salts was mostly oriented toward the bulk 
description (3). Particularly, the (2+l)D Chern-Simons 
term has been derived microscopically for the bulk effec- 
tive action of the system in Rcf . Q] . According to general 
principles (|||, a system where the Hall effect is quan- 
tized with an integer number N should have N gapless 
chiral electron states at the edges of the sample. How- 
ever, specific microscopic realization of such edge modes 
in Q1D organic conductors was not clear. In these sys- 
tems, the electron states at the edges parallel and perpen- 
dicular to the conducting chains must be very different 
because of the strong anisotropy. The edges parallel to 
the chains were studied microscopically in Rcf. 0. It 
was shown, that the +kp electrons on the first N chains 
at one edge and the — kp electrons on the last N chains 
at the other edge remain ungapped in the FISDW state. 
The velocity of these chiral modes is equal to the Fermi 
velocity vp along the chains. However, the edges per- 



pendicular to the chains were not discussed in Rcf. 0. 
Phcnomenological and numerical studies of multilayered 
quantum Hall systems did not address the micro- 
scopies and anisotropy of the edge states in Q1D case. 

In this letter, we develop a microscopic theory of the 
chiral states in Q1D conductors at the edges perpendic- 
ular to the chains. We find N branches of gapless edge 
excitations that are localized along the chains within the 
coherence length £jv = Tivf/ A^. Their group velocity 
perpendicular to the chains is v±_ = NA^b/h, where b 
in the interchain distance in the b direction. Thus, their 
localization length is much longer and the velocity much 
lower than the corresponding parameters Nb and vf for 
the edge states parallel to the chains. 



To explain the origin of the edge states |10|, it is use- 
ful first to consider a one-dimensional (ID) charge- or 
spin-density-wave (CDW/SDW) system. The mean-field 
Hamiltonian for the system can be written as 



n 



dxijr (x) 



2^2 



2m 



2Acos(2k F x + 9)^i'(x), (1) 



where we have chosen the x axis along the chain. ip{ x ) is 
the fermion (electron) field, m is the effective mass of the 
fermions, 2Acos(2kpx + 9) is the density-wave potential, 
and 9 is the phase specifying the density-wave position. 
To simplify the presentation, we suppress the spin struc- 
ture of the fermion field and the order parameter, which 
is not essential for our purposes. Following the standard 
procedure, we introduce a doublet of fermion fields, ip+ 
and "0-, with momenta close to the Fermi points ±kp: 



^(x) = ^ + { x ) e lkFX + ip-{x) e 



-ikpx 



(2) 



and linearize the energy dispersion relation near the 
Fermi points. The Hamiltonian can then be written in 
terms of the fermion doublet ip = , ip- ) as 



H 



dx^(x)(T Q e F - T z ihv F d x + r x Ae lT * e )ip(x), (3) 



where vf = kkp/m is the Fermi velocity, ep 
is the Fermi energy, and t x , t v , t z and tq are the 2x2 



h 2 k F /2m 



1 



Pauli and unit matrices acting on the fcrmion fields. 
From Hamiltonian (^), one obtains the eigenvalue equa- 
tion for the quasiparticlcs: 



-T z ihvpd x 



? ) i> = 



(4) 



where E is the quasiparticle energy as measured from the 
Fermi energy. For an infinite sample, Eq. (|4|) has well- 
known plane-wave solutions ip °c e%kx with the spectrum 



E, = ± 



Ti 2 vpk 2 + A 2 . 



(5) 



The spectrum (Q) forbids electron states within the en- 
ergy gap ±A. 

Now let us consider a semi-infinite CDW/SDW system 
occupying the positive semispace x > with an edge 
at x = 0. In this case, the wave function must vanish 
at the edge: i(?(x = 0) = 0, or equivalently ip+(x = 
0) = —i/j_(x = 0). With this boundary condition, in 
addition to the delocalized states (||), Eq. (||) also admits 
a localized state with the energy \E\ < A inside the gap: 



tp± = ±e 



-A cos ( 



k = -sin0/£. (6) 



The wave function (g) decays in the bulk at a length 
of the order of the coherence length £ = %v F /A. Eq. 
(^) is meaningful only when n > 0, thus the localized 
state exists only for tt < 9 < 2n. If the CDW/SDW is 
displaced by one period, the phase 9 changes from to 
2-7T, and the energy of the edge state (||) changes from 
A to —A. This process is the spectral flow |TTJ from 
the upper to the lower energy band (||). (The opposite 
spectral flow takes place at the other end of the sample.) 
The edge state @) is similar to the electron state localized 
at the kink soliton in an infinite CDW/SDW system jTJ||. 

Now let us consider a Q1D FISDW system. The Hamil- 
tonian of this system differs from Hamiltonian (|l|) by the 
interchain tunneling term with the amplitude |^Q: 

/dk / h 2 r) 2 

dx ip ^{x, k y ) ( - + 2t b cos(k y b - Gx) 

(7) 



2Acos(Q x x + 9) \i(j(x, k y ), 



where k y is the electron wave vector transverse to the 
chains. In the presence of an external magnetic field H, 
the interchain tunneling term contains the magnetic wave 
vector G = ebH/hc (c is the speed of light) because of 
the Peicrls-Onsager substitution k y — eA y /ch with the 
vector potential A y = Hx in the Landau gauge. Another 
difference between Hamiltonians (^) and ([?]) is that the 
FISDW wave vector Q x deviates from 2kp by an integer 
multiple N of G: Q x = 2k F - NG §§. To simplify pre- 
sentation, here we neglect the coupling t c between the (a- 
b) layers and model the system as a set of uncoupled 2D 
planes. Further, we only consider the electron tunneling 
between the nearest-neighboring chains and also assume 



that the transverse component of the FISDW wave vec- 
tor is zero. While these assumptions are too simplistic 
for computation of the quantities such as N, A, and T c , 
they do not affect our results about the edge states. 

Next we introduce the doublet of fermion fields ijj = 
{ip+, "0-) as m Eq. (|2|) and linearize the energy dispersion 
near ±kp. The Hamiltonian becomes 

/dk ( 
dx -p ip^(x, k y )\To€F - T z ihvFd x 

+ T 2t b cos(k y b - Gx) + T x Ae ir ^ NGx - e) ^ ip(x, k y ). (8) 

Then we make a chiral transformation of the doublet ip'- 

x/j(x,k y ) = e iT ^ x > k *^'(x ) k y ) ) (9) 

(10) 



(f(x, k y ) = (2tb/huj c ) sm(k y b — Gx), 



where Tioj c = TivpG is the cyclotron energy. In terms of 
the new field ijj' , the Hamiltonian can be written as 



dk 

dx ip''(x, k y )(T e F - T z ihv F d x 
Ztt 

+ T x Ae lT ' [NGx - e+2v(x ' k ^ ] )iP'{x, k y ). (11) 



The transverse hopping term is now transferred to the 
phase of the FISDW order parameter. Next, we expand 
that phase factor into a Fourier series 



J n (4:tb/huj c )e 



ir z n(k y b—Gx) 



(12) 



where J n is the Bcssel function of the order n. Substitut- 
ing Eq. ( |l2| ) into Eq. (p"l|), wc only retain the term with 
n = N in the series, since only this term does not have 
an oscillatory dependence on x and opens a gap at the 
Fermi level. This is the so-called single-gap approxima- 
tion |Tj|[l4|]. The Hamiltonian then becomes 

/dk / 
dx -p ip n (x, k y )^T Q e F - T z ihv F d x 

+ T x A N e^ Nk y b - e ^'(x,k y ), (13) 

where Ajy = AJpf(Ati,/hLo c ) is the modified gap ampli- 
tude. Comparing Eqs. ( |l3|) and (|3|), we see that the Q1D 
FISDW Hamiltonian has been mapped onto an effective 
ID CDW/SDW Hamiltonian with the gap amplitude Ajv 
and the order-parameter phase 9 ~ Nk y b. 

Now let us consider a semi-infinite FISDW system oc- 
cupying the positive semispace x > along the chains 
and infinite space in the y-dircction transverse to the 
chains. The wave function ip 1 must vanish at the edge: 
"01(0,?/) = —tl>'_(0,y). Following the same calculation as 
for the ID CDW/SDW, we find the edge states with the 
energies l^l < A^r inside the gap: 

V4(x, k y ) = ±e tk y y - KX /^, k = sin(Nk y b)/Z N , (14) 
E N {k y ) = -Ajy cos(Nk y b), £jv = Tiv f /An- (15) 
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FIG. 1. Energy dispersion of the electron states localized at 
the left (solid lines) and right (dashed lines) edges of the sam- 
ple as a function of the transverse momentum k y for N = 2. 



Here Epj(k y ) is the quasi-particle energy measured from 
the Fermi energy, and we have dropped the now unim- 
portant phase 9. The single bound state (jfy is now re- 
placed by the band ( |l5| ) of the edge states labeled by the 
wave vector k y perpendicular to the chains. These states 
( |l4| ) are bound along the direction of the chains, but are 
extended in the direction transverse to the chains. At 
the left edge of the sample, we require that k > and 
< Nkyb < 7T, and find A~ branches of the edge states 
in the transverse Brillouin zone < k y b < 2n. Com- 
plementary A~ branches of the edge states determined by 



the condition k < and it < Nk y b < 2ir exist at the 
right edge. The solid and dashed lines in Fig. [I] show the 
energy dispersions of the states localized at the left and 
right edges, respectively (for N = 2). It is clear that the 
group velocities of the edge modes dE^(ky)/hdk y have 
opposite signs for the left and right edges. Thus, they 
carry a surface current around the sample, as indicated 
by the arrows in Fig. ||. The sense of circulation is deter- 
mined by the sign of N, which is controlled by the sign 
of the magnetic field H . 

An effective action for the low-lying edge excitations 
can be obtained by linearizing the energy dispersion jl5| ) 
near the Fermi energy. Since there are N branches 
of these excitations, the effective action can be writ- 
ten in terms of an A^-component chiral fermion field 
4> = (fa, <fe, • • • , 0jv), as in Ref. (6): 



^edge 



dtdy$(t,y)(ifrd t + ihvj_d y )<f)(t,y), 
dE N (k y )\ _ NA N b 

' E N {k y )=0 " 



dk y 



(16) 



where t is time. The group velocity v± ( |16[ ) is same for all 
branches. The velocity is rather low, because it is propor- 
tional to the small FISDW gap Ajv- The activation en- 
ergy 2A N = 6K was measured in (TMTSF) 2 C10 4 |]§ at 
H = 25 T, where N = 1. The BCS relation A = 1.76 T c 
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FIG. 2. Sketch of the proposed time-of-flight experiment. 
The arrows indicate the directions of the edge modes velocities 
v± and vf- The thin lines indicate the conducting chains in 
the sample. The pulser sends a pulse, which is detected at 
different times t and t' by the detectors Dl and D2. 



gives even bigger value Ajy = 5.3 K for T c = 3 K at H = 
9 T in (TMTSF) 2 C10 4 @. Using the value A N = 3 K 
for N — 1, and the lattice spacing b = 0.77 nm, we find 
v±_ = 300 m/s. The group velocity for the edges perpen- 
dicular to the chains, v±, is three orders of magnitude 
lower than for the edges parallel to the chains, vf = 190 
km/s (see Fig. ||). Nevertheless, the total edge current 
/ is the same: I = (Ne/irh) j^ 2bN [8En / dk y ) dk y = 
NeAiy/irh = ev±/irb = 20 nA. For the typical sample 
dimensions L x = 2 mm and L y = L z = 0.2 mm and the 
interlayer spacing 1.35 nm, the edge current / produces 
the magnetic moment 1.2 nAm 2 and magnetization 15 
A/m. However, there arc additional contributions to the 
total magnetization coming from the edge states inside 
the energy gaps opened by the neglected terms in Eq. 
( |l2| ) below the Fermi level |l4|]. Magnetization has been 
calculated using the bulk free energy in Ref. |17| and 
measured experimentally in (TMTSF) 2 C10 4 m Ref. [@. 

The gapless edge excitations make a contribution C e 
to the specific heat per layer proportional to temperature 
T at T « A N : 
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Using the numbers quoted above, we find that t x = 
L x /vp = 10 ns and t y = L y /vj_ = 0.67 /is for N = 1. 
Thus, the dominant contribution to the specific heat 
comes from the edges perpendicular to the chains. Ac- 
cording to Eq. (|T7j), the specific heat must be discontinu- 
ous at the boundaries between the FISDW phases, where 
Ajv changes discontinuously |l7). Note that A^ cancels 
out in Eq. (|T7j). This result is in contrast with the phe- 
nomenological model of Ref. j^] , where the edge velocity 
was assumed to be the same for all FISDW phases and 
much stronger discontinuities C e /T ~ A^ were predicted. 
The bulk specific heat in the normal state is also 



3 



proportional to temperature: cj^ /T — TrL x L y /3hbvF- 
Its ratio to the edge specific heat (|l^) is roughly equal 
to the ratio of the volumes occupied by the bulk and 
edge states: C { b n) /C e = L X /2£ N = 2.1 x 10 3 . In 
the FISDW phase, the bulk specific heat C ( b F) /T « 
^18/^ (Aw/? 1 ) 5 / 2 exp(-AAr/T) cf ] /T is smaller than 
the edge one at the sufficiently low temperature T < 
Aiv/14. Such a regime could have possibly been achieved 
at H = 9 T in the experimental measurements of 
the specific heat |l6| performed at T = 0.32 K in 
(TMTSF) 2 C10 4 . We estimate the absolute values of 
C e /T as 2.2 x 1CT 3 mJ/(K 2 mole) and C ( b n) /T as 4.4 
mJ/(K 2 mole). The latter estimate roughly matches the 
experimental value 5 mJ/(K 2 mole) Jl(J. 

The edge states picture was confirmed in a multilay- 
ered GaAs system by demonstrating that the conduc- 
tance G zz perpendicular to the layers is proportional to 
the number of the edge states N and the perimeter of 
the sample ]l9| . A similar experiment on the Bcchgaard 
salts produced inconclusive results pO}j . As a definitive 
proof of the edge states picture for a FISDW, we propose 
a time-of-flight experiment analogous to those performed 
on GaAs in Rcf . j2^] . The experimental setup is sketched 
in Fig. ||. A pulse is applied by the pulser at the center of 
the edge perpendicular to the chains. The pulse can be 
electric, as in Rcf. pi] ], and produce a charge perturba- 
tion in the occupation number of the edge states. Alter- 
natively, the pulse can be thermal and perturb the Fermi 
distribution function of the edge electrons, or magnetic 
and change the spin population of the edge states [f22| . 
The pulse will travel with the slow velocity v± along the 
edges perpendicular to the chains and the fast velocity 
vp parallel to the chains. The pulse will reach the detec- 
tor Dl at a time t = L y /2v± and other detector D2 at a 
time t' = 3t + 2L x /vp sa 3t. The difference between t and 
t' is a signature of the broken time-reversal symmetry in 
the system. The inverse delay time 1/t is proportional 
to v±, given by Eq. (|l^), and has discontinuities at the 
FISDW phase boundaries due to discontinuity of both N 
and Ajv- The pulse must be shorter than t = 0.33 /is for 
clear resolution and longer than %/ A at = 2.6 ps, so that 
only the low-lying excitations are probed. 

In conclusion, we have shown, starting from a mi- 
croscopic Hamiltonian, that there are N chiral gapless 
electron states bound to the edges perpendicular to the 
chains in a FISDW phase characterized by an integer 
N . The gapless modes contribute a term linear in tem- 
perature, which dominates the specific heat at low tem- 
peratures. The group velocity of these excitations is 
rather low: v± = NA^b/fi as 300 m/s. We have pro- 
posed a time-of-flight experiment to probe the propagat- 
ing edge modes directly. The chiral character of these 
states makes them similar to the edge states in super- 
conductors with broken time-reversal symmetry, such as 
d + is and d x 2_ y 2 + id xy singlet superconductors and 



Px + ipy triplet superconductors Si'2Ru04 and 2D 3 Hc-A 
|||. (TMTSF) 2 X in the superconducting state is ex- 
pected to have nonchiral midgap edge states [E4j . 
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